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Abstract 

We demonstrate how to obtain explicitly the propagators for quan- 
tum fields residing in curved space-time using the heat kernel for which 
a new construction procedure exists. Propagators are determined for 
the case of Rindler, Friedman-Robertson- Walker, Schwarzschild and 
general conformally flat metrics, both for scalar, Dirac and Yang-Mills 
fields. The calculations are based on an improved formula for the heat 
kernel in a general curved space. All the calculations are done in d = 4 
dimensions for concreteness, but are easily generalizable to arbitrary 
d. The new method advocated here does not assume that the fields 
are massive, nor is it based on an aymptotic expansion as such. 
Whenever possible, the results are compared to that of other authors. 



1 Introduction 

The calculation of propagators of various fields in a curved background is of 
utmost importance in theoretical physics, since it essentially provides the only 
way of calculating scattering processes in, say, the early universe or around 
a black hole. As there is at present no fully satisfactory theory of quantum 
gravity, the only way of studying quantum effects in the presence of strong 
gravitational fields is through quantum field theory in curved spacetime, 
what one might call semi-classical gravity. The quantities of interest there 
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are first and foremost the propagators, the effective actions and the energy- 
momentum tensor. Therefore a lot of work have been done in the past two 
or three decades in this field, see e.g. [I], 0, [3|, [|]. In this paper we will 
concentrate on the calculation of the propagators of quantum fields of spin 
zero, one-half and one. To the best of our knowledge this is the first time 
propagators of non-zero spin fields is determined in a general background. We 
will leave the study of gravitons on arbitrary backgrounds for future research 
- the main aim of this paper is quantum field theory in curved spacetime. 
In previous papers f|, || we have developed a method for the determination 
of generating functionals of quantum fields residing in curved space-time. 
This was done using the heat kernel method and fairly reliable methods for 
determining this quantity were devised 0. One could decide to determine 
propagators for the quantum fields from the generating functional but it 
turns out that one can construct the propagators directly from the heat 
kernel and this latter approach is the one pursued in this paper. It should 
be emphasized that the resulting propagator is not just the free or bare 
one, as the calculations put forward here are indeed non-perturbative, at 
least in the coupling to the background fields (usually just the gravitational 
field, but occassional comments are made on how to include other kinds of 
backgrounds, e.g. external Yang-Mills or Higgs fields). 

In section 2 we relate the heat kernel of a differential operator A, associated 
with the quantum field in question, to the corresponding propagator. Then 
we show how to explicitly determine propagators using previous results on 
constructing the heat kernel. This is just the Schwinger-DeWitt proper time 
formalism, but we use an improved expression for the heat kernel. Section 3 
is devoted to a discussion of the corresponding vacuum state and a discussion 
of the Hadamard condition. We exemplify the approach in section 4 where 
we also compare our results to that of other authors and finally provide a 
brief conclusion and outlook in section 5. 
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2 Determining the Propagator From the Heat 
Kernel 

Call the operator of interest, which generally varies in space, A = A(x), and 
the corresponding eigenf unctions ip\ so that 

Atfj x = A^a (1) 
The heat kernel Ga{x, x', a) is the function satisfying the heat kernel equation 

d 

AGa(x, x ; a) = ——Ga(x, x ; er) 
da 

subject to the boundary condition Ga(x,x';0) = 5(x,x'). 

Note that this equation is satisfied by (the spectral representation of the heat 

kernel) 

G A (x, x', a) = J2 M^Wx(x')e^ (2) 

A 

The Green's function on the other hand is a solution to the equation 

AG(x,x') = 5(x,x') (3) 
which is satisfied by (the spectral representation of the Green's function) 

GforO^E^ACO^OOA"" 1 ( 4 ) 

A 

Put together with equation (0) this yields the relationship between the 
Green's function and the heat kernel: 

G(x,x') = — / cI(jGa(x, x 1 ] a) (5) 
J o 

(provided that A 7^ 0). This equation alone shows that the usual asymptotic 
expansion, due to Schwinger and DeWitt, which is only valid for a — > 0, can- 
not be relied upon to provide good propagators. Furthermore, for massless 
fields, the Schwinger-DeWitt expansion would lead to a manifestly divergent 
cr-integral. In the following section we will write down an improved expan- 
sion for the heat kernel valid for o large as well, and which moreover leads 
to closed expressions for the expansion coefficients. This improved expan- 
sion is convergent even for m — 0, and the method of calculation is readily 
generalized to higher spin too. 
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2.1 Heat Kernel and Propagators for a Scalar Field 



The operator involved in the case of scalar particles is the (curved space) 
d'Alembertian (minimally coupled scalar fields) or the d'Alembertian plus 
some function (non-minimally coupled scalar fields). So we express the 
d'Alembertian in terms of the local vierbeins, ef, = l 2 ^- where Greek indices 

' M ox^ 

refer to general coordinates while Latin indices refer to the local inertial 
frame (the metric is g^ u = rjabe^e^ where r] a b is the Minkowski metric, g is 
the determinant of the metric and e the vierbein determinant, e = y/g)' 

= -d,(e V ab e^ b d u ) 
e 

= -e^d m (ev ab e^ b e n M 
e p 

= h™d m (ev ab e^ b e n u )d n + r) ab e™e^ b e n „d m d n 

= □o + ^(<9 m ( e <))3 a (6) 

where D = r] ab d a db is the d'Alembertian of a comoving observer the heat 
kernel of which is known to be @ 

A.(x x' ) 

G (x,x' ; a) = (4ira)~ 2 e (7) 

where A(x,x f ) = (J^ds) 2 , is the geodesic distance squared (i.e. half the so- 
called Synge world function) which is just (x — x') 2 in Cartesian coordinates.^ 
We will refer to O as the "flat" d'Alembertian, and consequently to Go as the 
"flat" heat kernel. Proceed to remove the first order term of the Lagrangian 
by the substitution 

Go = G(x, x'; <r)e-3 / (g) 

1 The heat kernel (u\) is the straightforward covariant generalization of the result in 
Minkowski space-time G = (Aito)~ 2 cxp(— ^ )■ In d dimensions one simply has to 
replace (4na)~ 2 by (4iT<7)~ d / 2 . The vierbeins are introduced to make this transition from 
flat to curved space more transparent. 
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The integral in this expression is always easy to calculate as the term dif- 
ferentiated is just the reciprocal of the term with which it's multiplied by 
so that one gets a logarithm of a product of vierbein components which in 
turn makes the exponential the reciprocal of the squareroot of this product. 
It turns out, however, that we actually do not need the explicit form of this 
integral as it cancels out in the final expression. 

When making the substitution (||) in equation (H) one gets a zero'th-order 
term which must be added to the one that was there in the first place (i.e. 
to £-R). One thus gets the following heat kernel equation (£ = in the case 
of minimal coupling): 



1 /l _ \ 2 1 (\ 



G(x, x'; a) = —d a G(x, x'; a) 
(9) 



written compactly as 

(□ + f (x))G(x, x'; a) = -d a G(x, x'; a) (10) 

with 

fo = ZR+ \ (e-%{ee£)) 2 - V {e~%{e^)) (11) 

The method for solving this equation devised in papers makes use of the 
following trickF] 

G(x,x';a) =G (x,x';a)e- T ^ x '' a) (12) 

Taylor expanding T 

oo 

T(x,x';a) = £ r n {x, x')a n (13) 

n=0 

one gets the following recursion relation for the coefficients 

n— 1 j 

nr n = -n r n -i + dr n' ■ di~ n -x-n> + -dA ■ dr n (14) 

n'=0 1 



2 One should note that this expansion is not the usual asymptotic one due to Schwinger 
and DeWitt (see e.g. 0]), but actually holds for all values of a. Furthermore, as the 
recursion relation will show, it is actually rather easy to evaluate the coefficients, especially 
along the diagonal x = x'. The expansion shown here works for massless fields aswell, 
contrary to the Schwinger-DeWitt expansion. 
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with To = +| / e~ 1 d fl (eefi)dx n , as follows from the boundary condition 
linio-^o G(x, x'; a) = 5(x,x'). Along the diagonal x = x', this recursion re- 
lation simplifies and the result is listed in J5J. In the general case x 7^ x', 
however, the coefficients develop a dependency upon (the values of the cur- 
vature and its derivatives along) the geodesic from x to x' (assuming the 
existence of such a curve - otherwise A has to be put equal to 00, leading to 
zero amplitude for the propagation). Explicitly, 

r n {x, x') = e -2rV*(*.*') f 1 ^ n _ 1 (r)e 2nr dr + c n e~ 2n ^ ^ x ' x ' ] (15) 



where the integral is along the abovementioned geodesic, and the functions 
# n _i are given by 

n 

g n (x,x') = -a r n (x,x') + Yj dr n> ■ dr n _ n/ 



n'=\ 



with the first few functions being 

So = h = f + d n (e~%e^+e- 2 V mn d,(ee»Jd u (ee:) 

91 = - D oi~i(x, x') + 28t ■ dri 

92 = -a r 2 (x,x') + (dT 1 ) 2 + 2dT -dr 2 



and so on, with f given by (|TT|), g is not just simply f , as one might 
otherwise expect, a fact due to the boundary condition G(x, x'\ a) — ► 5(x, x') 
for cr —>■ 0. The variable r is the proper-time of a comoving observer, scaled 
such that at r = the observer is at position x, whereas at the later time 
t = 1 (s)he is at x'. 

Inside the integral along the geodesic the differential operators d, O simplify 
to become 



d d 
d . d 2 



d P = ^ = e p %(r)- = x p (r)- (16) 



□0 = V Pq d P d q = r ] ™ Xp x g — + x 2 — (17) 

where x^{r) describes the geodesic. These differentiations are consequently 
rather straightforward to carry out. 
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The coefficients c n are independent of the curve, they only depend upon the 
value at the end points, in fact 

lim T n (x, x') = c n (x) and lim r n (x,x') = c n {x') (18) 

x'—>x x—*x' 

and they are thus equal to the coefficients found in H. They satisfy 



n 



+ l)c n+ i = -n c n + dc n' ■ dcn- 



(19) 



n'=0 



and thus the relevant ones aref] 

ci = f 
c 2 = 



c 3 = inj»/ o _i(0/ o )* 
- -i(5/o) 2 

where we have decided to only include first and second derivatives of /o 
(essentially the curvature) in the last approximation. These expressions are 
for massless fields, masses are accomodated by adding m 2 to f , and therefore 
to Ti,c\. Written out in full, the first non-trivial coefficient function, g\, is 
then 



9i 



A- 3 / 2 (9A) 2 + A-^idA) 2 - A~V 2 n A 



Tl + 



, p Ci 
+ 



-2\/A 



A-^(^A) [jf d p {g e 2T )dr + -d, 

Jo 2 

V mn e-%(ee») [-A~ 1 / 2 (a m A)r 1 + e~ 2VE d m (foe 2T )dr + d m c x 



where all derivatives outside the integrals are with respect to x and not x', the 
derivatives inside the integrals, however, are to be understood as in ( p~6| - p~7| ) . 



3 Strictly speaking the coefficients are \{c n {x) + c n (x')), but for simplicity we have 
decided not to write them in this way, thus a symmetrization in x, x' is always understood. 
In the particular examples this is carried out explicitly. 
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The next coefficient function, g 2 , becomes 



92 



A-V2( 5A )2 + 4A- 1 / 2 (9A) 2 - 2A- 1 / 2 D A 



t 2 + 



4A- 1 / 2 (9 p A) [jf d p { gi e AT )dr + -d p c 2 
Jo 2 

2A- 1 / 2 r 1 e- 2 ^(9 p A) d p (g e 2T )dr + ^ pCl ] } 



+ 



-2A- 1 / 2 (9 m A)r 2 + e~ 4 ^ { jf* 9 m (^e 4r )rfr + <9 m c 2 



As is apparent, the expressions quickly get somewhat involved. Note, how- 
ever, that most of the needed operations are trivial - most of the terms are 
simple derivatives of A and fo- The only non-trivial part being the integral 
along the geodesic. By the same token, we will usually only list fo and t\ in 
the examples to follow. 

One should notice that in principle one could find the coefficients to any 
chosen order, we will however stick to the approximate solution with r n « 
, n > 4, corresponding to d 3 R, d 4 R, ... being negligible. In fact we can 
make do with the even more tractable form (which is adequate^) 

G(x,x';a) « G (x, x'; a)e- T1CT {l + r 2 {x, x')a 2 + r 3 (x, x')a 3 ) 



(Ana 



-2, 



X,X I I l\ 

-L_L- T1 ( X ,X W]_ _|_ 



T 2 {X 1 X )<J' 



+ r 3 (x,x , )a 3 )(20) 



where expression ([?]) for the "flat" space heat kernel has been inserted. 
Other modifications of the Schwinger-DeWitt asymptotic expansion exists. 
In fact Parker and coworkers, @, have suggested that a factor exp(— £-Rcx) 

4 As general relativity is only renormalizable upto one loop order, it would not be mean- 
ingful to go to higher order in this approximation; the higher order terms are essentially 
third and higher order derivatives and powers of the curvature, corresponding to a higher 
loop order. This approximation corresponds to a regime in which R ^> (dR) 2 , O R, i.e. 
strong but slowly varying curvature. The other extreme (R <C dR, i.e. strongly varying 
curvatures, e.g. in the space-time foam) can be treated similarly, the coefficients c n will 
then be c„ » (-l) 71 ^ 1 ^a^ 1 f Q , and similarly for g n ,T n . 
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was present if one performed a partial summation of the usual expansion. 
The basic new thing in our approach is the actual possibility of a systematic 
calculation of all the coefficients in the expansion - in fact we've got a gen- 
eral, closed expression, (|1^). The expansion of Parker et al. differs from ours 
in that whereas we expand T in exp(T), they expand exp(T), thus getting 
simply modifications of the Schwinger-DeWitt coeffcients. The main point 
of the approach put forward in this paper is the actual abbility to systemati- 
cally find the expansion coefficients both for x 7^ x' (the general case) and for 
x = x' (which turns out to be very simple). It is the use of the exponential 
which makes this simplification possible. 

As pointed out by Ford and Toms, ||, there is a genuine need for non-local 
terms in the heat kernel (and hence a need to go beyond the Schwinger- 
DeWitt expansion) and as it turns out, the coefficients obtained by the 
method put forward here are in fact non-local. Their geometric meaning 
is even rather transparent as they involve the integration along a geodesic of 
various curvature-related quantities. 

The heat kernel can also be written as a path- integral, |9[], and the usual Van 
Vleck-Morette determinant occurs when one considers geodesies, [[H| . 



The cr-integrals involved in determining the Green's function from equation 
(5) are now of standard type leading to the following result [[11 



G(x,x')^ (47T)~ 2 2 . 7 /J K_ 1 (y/A(X,X')T 1 (X,X')) + 



T 2 {X,X) 




v 4ti(x, x') 



K 1 (^A{x,x')t 1 (x,x')) - 



T 3 {X,X 



/)^^K 2 (VA(x,*0ti(x,*0)) (21) 



where K n is a modified Bessel (or MacDonald) function which has the fol- 
lowing series expansion (needed in the case of gauge bosons treated below 
and for the proof of the Hadamard condition being satisfied): 

«■.(») = (-d" +1 t ('» r - 5*<» +k+1) ~ + 1) ) (f 

n-l /„ _ u _ i\| / rr.\2k-n 
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with ip(x) = -£-lnr(x) is the Euler psi (or digamma) function. 
Using the recursion relation | I 

zK n _i(z) - zK n+1 (z) = -2nK n (z) 

together with K_ n = K n , one can reexpress ( f21~D in terms of Kq and K\ 
-±K , 0, only as 



G(x, x') 



-a (x,x')-^- + (3 (x,x') 
dz 



K (z) = a (x, x')Kx (z)+p (x, x')K (z) 

(23) 

with z = J A(x, x')ti(x, x'). To the chosen order we have 

a (x,x') = -2(47r)- 2 ( ,/^ + r 2 ,/A + T3 /" ) {21) 




Po(x,x') = --M47T)- 2 (25) 



Note that the form (|23[) of the Green's function is exact, it is only the explicit 
expressions ( p4"| - f2"5| ) for «o,/5o which are approximate. 

Thus the problem of finding the Green's function of a scalar field in some 
geometry, is reduced to evaluating the two functions fo, fo (and from these 
r n ), containing the information about the metric structure. The calculation 
of these is, as will also be apparent from the examples to follow, rather 
st r aight f or war d . 

We will see that also the propagators for fields with non-zero spin can be 
written in this form, involving only the two modified Bessel functions Kq, K\. 
The coefficients will contain all the spin-specific data, and will be denoted by 
a s , (3 S for spin s. Only for s = are these scalar functions, for bosons of spin 
s they will in general be 2s-tensors, whereas for spinors they will take values 
in the Clifford algebra (for Rarita-Schwinger fields, s = 3/2, they would carry 
two Lorentz indices besides their two Clifford algebra (i.e. spinor) ones, and 
so on). 



2.2 Heat Kernel and Propagators for Spin 1 Gauge 
Bosons 

The heat kernel for spin 1 gauge bosons has been determined in || and 
because it is the result of a lengthy calculation we simply quote the result. 
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The heat kernel equation can, in the mean field approximation chosen in 
be written asR 



8 a h 8™d p dP + 5 a b (d n e m » - d m e») e p d p + {gf b a c (d n A mc - d m A c J) 



+ (-5™g 2 f ebc frA e p Av d ) 



Gt(x,x':a) 



6b da 



(26) 



where G = G q m is a matrix valued function. Defining 

£ rnp = ( dne m» _ Qm ^ g p ^ 

one can remove the first order term in equation (|26"D by substituting (the 
quantities in the below equation are matrices) 

G = Ge--* E (28) 

where d p E™ = £™ p . The heat equation then becomes 



5™d p d p G q m — 



\d v ET + ~C^4 + (gf b a MA mc - d m A c n ))+ 



(^n9 2 fe b cf d a c A e p A pd ) 



d ~ 
ni = —Hi 



(29) 



an equation of the same functional form as for the non-minimally coupled 
scalar field. The solution is then (in matrix notation) 

G{x, x'; a) = Gu (x, x'- a^-- 4 * 7 ^ 2 -^ 3 



withg 



B" 



C 



-2V~A 



A™(x(T))e 2T dT + A™e~ 2V7£ 

2 jo 

1 -gVa [\d AT)(d p A k n )e 6r dT + C™e~ 6 ^ 

3 Jo 



(30) 

(31) 
(32) 
(33) 



5 In || we have calculated the effective action using the heat kernel without using a 
mean field approximation, but the result from which we use here, is conceptually 
simpler. 

6 Even though the matrices A, £>, C do not in general commute, there is no problem with 
this solution, as their commutator will inevitably be a higher order term. But the formula 
is not as useful as for the scalar case due to the complications of the gauge field. On the 
other hand, the couplings and mean field values need not be small, as the result presented 
here is non-perturbative. 
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where 

j% = d p Sr + \srS k np +{gf b a c {d n A mc -d m A c n )) + 
(^g 2 f ebc f d a c A;A^ d ) 

K = 

corresponding to Ci = f , c 2 , c 3 respectively, and 

□o = <9 P ,9 P = = grd^dv + rT^^R (37) 

as before, the heat kernel of which is given by equation ([?[). Analogously to 
the scalar case, the quantity A is essentially a curvature term. The approxi- 
mation involved in the evaluation of C, C above is, as for the scalar case, the 
omission of a fourth order derivative in the curvature, namely O^A. 
A few comments concerning the actual calculation of this quantity are in or- 
der. Now, if the matrices A, B, C were simply diagonal, one could apply the 
Bessel functions directly to the diagonal values, which would then become 
scalar boson propagators. Even when A is diagonal (which holds for abelian 
fields), it will not a priori commute with B, C and these can therefore not in 
general be diagonalized, leading to some technical difficulties in the practi- 
cal application of this formula. In practical calculations it will of course be 
tempting, and, for most purposes, good enough, to ignore the B and C terms 
of the heat kernel (which is strictly speaking only good enough for abelian 
or weak fields and couplings) or to resort to numerical calculation of these 
higher order terms. This is unsatisfactory (and, moreover, rather cumber- 
some), but we have not yet been able to come up with a better solution. 
Referring back to equation (||) the propagator becomes (compare with the 
scalar case of equation^!])) 

D™(x,x')= - (4tt)- 2 (2(^)^/2(^/2^ (K^(VAA)) h n + 

l -C% -A^T 1 )* (K 2 (VAA))) l n ) (38) 
12 



(34) 

(35) 
(36) 



which we will write as (in matrix notation) 



D(x, x') = ai (x, ar^iCV^n) + x')K (^A^) (39) 
To the chosen order 

a* = 2(47r)- 2 ((^A)- 1 / 2 A 1/2 + (^A) 1 / 2 BA- 1/2 + ^ACA- 1 ^ (40) 

/3i = ^AiAny'CA- 1 (41) 

with the Lorentz indices suppressed. 

In order to determine the Green's function for a Yang-Mills field, all we have 
to do is to evaluate the matrix A, and from that the matrices B, C. Just like 
the functions f , f for the scalar case, these matrices only depend upon the 
vierbein (and the mean field, which can be expressed in terms of the vierbein 
as in H), but unlike the spin zero case, their evaluation is somewhat more 
complicated (but still feasible). 
Lastly we have to consider the case of fermions: 



3 Heat Kernel and Propagators for Dirac 
Fermions 

Consider the heat equation for a Dirac field in curved space-time 

d 

(?V - m)Gi/ 2 (x, x; a) = -—Gy 2 (x, x; a) (42) 



where[] 



= e^ a (d, + uy bc ) (43) 



7 It would be very easy to include, say, a coupling to a Yang-Mills field igA^^Tk, 
e.g. using the mean-field derived in the previous section, or a Yukawa coupling to a scalar 
or pseudo scalar field. The latter two would just amount to the addition of a term to 
the mass, whereas the Yang-Mills coupling would appear in an extra contribution to the 
function F a to be introduced below. The Yang-Mills field, however, would lead to a non- 
diagonal form ofy 2 , as it would now contain a term like a^ v F^ v Tk, i.e. a spin-magnetic 
field coupling. The introduction of a background torsion would lead to a similar term as 
shown in S . 
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This equation can be rewritten as 

{i^ida + F a ) + i l5l a F a - m)G 1/2 = (44) 

where 

F a = 8te»u b ;5 d b Vac (45) 
F a = -Aie^s bc d a (46) 

Here we have used that the a^ U) the generator of SO(l, 3)-transformations 
when acting upon Dirac spinors, is an element of the Clifford algebra, a mn = 
\i [lm->ln[- The spin connection can be written in terms of derivatives 
of the vierbeins. The heat kernel for a Dirac spinor takes values in the 
Clifford algebra, and can thus be expanded on the sixteen dimensional basis 
1) 75) 7a ; l5la, a ah- Due to this complication we have not been able to find 
an expression for the heat kernel for a Dirac fermion. We can still find the 
propagator, however, by noting 

S(x, x') = (ty + m)G'(x, x') (47) 

is an inverse of the Dirac operator provided G' satisfies 

{tf - m){tf + m)G' = (-y 2 -m 2 )G" = 1 (48) 

Noting that 

Y 2 = (□ + 0/2)1 (49) 
where 1 is the unit element in the Clifford algebra (i.e. in this particular 
case a 4 x 4 unit matrix) we see that G 1 is the Green's function for a non- 
minimally coupled scalar field. Had we considered a spin-3/2 fermion (a 
Rarita-Schwinger field) we would have made the same Ansatz but with G' 
replaced by D', a similar propagator for a spin-1 boson. 
The propagator then becomes 

S(x,x') = + m){a Q K 1 +(3 Q K Q ) 

= a xj2 K x + p 1/2 K (50) 

where cxi/2, P1/2 are Clifford algebra-valued. To the chosen approximation, 
these new coefficient functions are 

Oii/2 = ffl + m)a - i(z~ 1 a + Pojfz (51) 
A/2 = («y + m)/3o - ia$z (52) 
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with 

z = y'Ari (53) 

We have now written S(x, x 1 ), the propagator for a Dirac field, on the same 
form as that of G, the scalar field propagator. 

3.1 On the Corresponding Vacuum and the Hadamard 
Condition 

Now, in order for a two point function to be a propagator, it must be possible 
for it to be written as the expectation value of the (path ordered) product of 
two field operators with respect to some vacua (0 ou t|, |0m)- As is well-known, 
[0, |2|, |3|, the concept of a vacuum is somewhat complicated in a general 
curved space-time background, and different definitions will in general lead to 
different results. Some particularly important vacua are the asymptotic ones 
(assuming the space-time manifold to be asymptotically flat), the adiabatic 
(assuming sufficiently small or slowly varying curvature) and the conformal 
one. Other vacua, relying on the analogy with thermal field theory, 0, |], 
are the Unruh and Hartle-Hawking ones, the applicability of these is usually 
restricted to Rindler and Schwarzschild space-times, however. It turns out, 
by construction, that our definition of a Green's function entails a natural 
vacuum. Recall the spectral representation of the Green's function 

A A 

In a second quantized formalism, these eigenfunctions are of course to be 
replaced by operators, ip\(x) = ip\(x)\O x ). Remembering that we're relying 
on comoving coordinates when determining G, we see that our vacuum is 
that of a freely falling observer, and is thus well-defined on any manifold, so 
there is no need to impose global constraints on its topology. The propagator, 
however, will vanish whenever x and x' are not causally related. Due to this 
relationship with freely falling coordinates we will refer to our vacuum as the 
freely falling one, or, somewhat more amusingly, the "elevator vacuum." 
The question of when a candidate for a two point function is actually a 
physically valid propagator, has been studied intensively by Bernard S. Kay 
and coworkers, see |3J for a review. They impose what they call the Hadamard 
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condition, which, in the heuristic formulation of DeWitt and Brehme |12[] (see 
also Wald |J), states that for X ~ X 

JA vvm (x,x') ( 4 1. ,, . ,.\ 

G{x, x ) ~ J — — + u(ar, x ) In -A(x, x ) + w(x, x ) 

67T Z \A{X,X') 2 J 

(54) 

with Awm the Van Vleck-Morette determinant 

A VVM {x,x') = -(g(x)g(x , ))- 1 / 2 det(~W IM V u A(x,x')) 

and where v, w are smooth functions as x' — > x. Kay and Wald, Q, show that 
this is a physically reasonable requirement, and give a much more rigorous 
formulation of it. We want to show that our candidate, G(x,x') (and hence 
also D^(x, x') and S(x,x') for spin ^ 0), meets this condition. 
Now, writing G(x, x') in the compact form 

G(x, x ) = -2(Att)- 2 [a(x, x^K^z) + (3(x, x')K (z)\ 

with z = y / A(x, x')ti(x, x') and with a, /3 as given earlier in terms of A, r n , 
and using 

K (z) = -ln^ + V(l) + 0(z) 
^(z) = 2^ 1 + 0(1) 

we get 

G(x, x') = a(x, x) A -1 — /?(x, x') In - A + 7(x, x') (55) 

where 

- 2(47r)" 2 a(x, x') = a(x, x')A~ 1/2 + finite terms (56) 

and 7(x,x') is well-behaved as x — > x' as is a(x, x') and /3(x, x'). Our can- 
didate thus has the right form, and hence meets the Hadamard condition. 
Notice, however, that we escape the logarithmic divergence as (3 = 0(A), 
i.e. while we do have a logarithmic term it does not lead to a divergence as 
x' — > x. This is a most fortunate characteristic of our vacuum. 
The way one usually finds a propagator, |]J, is by using the positive frequency 
solutions to the equation of motion - the restrictions on the frequency being 
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there to ensure causal propagation. With the method put forward here, this 
requirement can be relaxed as we no longer need to know a complete set of 
solutions. Causality, moreover, is taken care of not by a process of time or- 
dering, which might be difficult to define rigorously in a general background 
over large distances, but in terms of path ordering, or, if one wishes, in terms 
of the proper time of the comoving observer. We thus get a causal propagator 
by restricting attention to time-like (for massive fields) or null geodesies (for 
massless ones). 

We will make a few comments on the comparison of our propagator with 
that of other authors, if such result exist (to the best of our knowledge), 
and it will turn out, that the results presented here disagree with all other 
calculations, due to the different notions of vacua. As our propagator does 
satisfy the Hadamard condition it is physically valid. Moreover, it is more 
general and apparently higher-order. It has for instance been argued that an 
adiabatic vacuum for spatially flat Friedman- Robert son- Walker geometries 
is only a Hadamard state if it is taken to infinite adiabatic order, |24]. 
Let us finally make a comment on the Schwinger-DeWitt expansion, 



oo 

A 



G(x, x'] a) ~ (4vra)- 2 e- m a n(%, s>" 

n=0 

Our expansion can be written in a similar form, by simply Taylor expanding 
the exponential. Denoting the coefficients by b n , we have 

b = 1 
b\ = ti 

b 2 = -Ti + T 2 

and so on, but b n 7^ a n in general, a difference stemming from the different na- 
ture of the two expansions. At a first glance, our coefficients seem to depend 
only on the curvature scalar and not on the Ricci or Riemann-Christoffel 
tensors, but actually these tensors appear through the derivatives of the 
geodesic distance squared. It turns out, as we have emphasized earlier, that 
it is exactly this use of an exponential that makes a closed expression for the 
expansion coefficients possible. An expression, moreover, valid for the gen- 
eral 7^ x' and reducing to a rather simple recursion relation for x = x'. 
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This latter case is important for practical calculations of, say, the renormal- 
ized energy-momentum tensor, [|], ||), and with the usual Schwinger-DeWitt 
expansion or with the modification due to Parker et al, only the first few 
can be found. On the contrary with the use of the exponential expansion 
suggested here, one can in principle (and even in practice) systematically find 
the coefficients to whatever order one wishes. We have used this possibility 
to calculate effective actions and energy- monentum tensors elsewhere || , and 
further work is in progress on the applications of this expansion. 
When studying coincidence limits, it is useful to evaluate A for x, x' infinites- 
maly close. We thus concider two points x, x' with x' = x^ + et M (in practice 
only one component of t M will be non-zero). We can then approximate the 
geodesic by a straight line with a small correction. Letting, for a non-null 
geodesic, £ = t^, as is common in the point splitting method, [|, we get 



T"2 



r 3 



VX)fo(x) + ~ePdJ {x) 



-(e 2 + l)f (x) + e 
^°o/o(aO(l " WX) + ^(e 4 - l)et»n d,nMx) - 

\^d,uMx) + o{t 2 ) 



+ 0(e^7) 



(58) 



144 



hl-6eVz)(dJ (x)) 2 + O(e 2 ) 



(59) 



A quick application of the propagators found in this paper is the cal- 
culation of mean fields. As we have just proven G(x, x') is divergent as 
x' — > x, but we can use (^-function regularization to obtain a finite value 
for < <f)(x) 2 >= G(x,x), the mean field. Going back to the expression of 
the propagator as an integral of the heat kernel, equation (|5|), then we see 
that the x 1 — > x limits gives a T(— l)-singularity. Such divergencies are to be 
removed by taking principal values, [25fl. We therefore get 



(0(x) 2 ) rcg = (4tt)- 2 -( 7 - l)n + t 2 t^ + r 3 rf 2 + ... 



(60) 



as the final result. 
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4 Examples 



To illustrate the technique developed above, we have chosen a few important 
examples. The first is Rindler space-time, possibly the simplest non-trivial 
example, and used for instance as an approximation to the Schwarzschild 
solution. The second and third examples are probably the most important 
ones from a purely theoretical point of view, as they are the classical solutions 
of Friedman- Robert son- Walker and Schwarzschild respectively. We finish 
off with a general conformally flat space-time (e.g. de Sitter space). For 
simplicity, we will just give the case of minimally coupled, massless scalar 
fields. Furthermore, we will just list the simplest expressions, i.e. t±,A and 
the coefficients Ci, c 2 , c 3 , A, B, C, the remaining terms (including all those for 
the fermions) follow by suitable differentiations, which, although very simple 
to carry out, leads to rather involved expressions. 



4.1 Rindler Space-Time 

The Rindler space-time is given by the line element 

ds 2 = (gz) 2 dt 2 — dx 2 — dy 2 — dz 2 



(61) 



where g is some constant, it is conformal to a wedge in Minkowski space-time 
and represents an accelerated observer in flat space. With this the vierbeins 
can be chosen to be 



( gz\ 








f ^ 




(°\ 








1 



e a — 



1 


e a — 






V o J 




V o J 




V o J 




V i / 



(62) 



A geodesic in this space-time takes the form 

t'-t 



t(r) 
x{t) 

v(t) 

z(r) 



sinha 2 

(x' — x)t + x 

(y' -y)r + y 



sinh — + t 

a 1 



aw In cosh — + b 



(63) 

(64) 
(65) 

(66) 
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with a, b two constants given by 

aVb = z a-s/lncosha -2 + b = z' (67) 

Thereby 

-i r _ / t' -t \ 2 

g 2 a 2 (In cosh ra 2 + 6) ( — cosh 2 ra 2 — (x — x) 2 — (y' 

\ smh a~ z I 



o 



It t 
•- tanh — -(In cosh — + b) 
4 a~ 2 a 2 



-1/2 



1/2 



which we have not been able to reduce any further. In practice, though, 
since the integration is over a compact interval, it is fairly easy to carry out 
numerically 

4.1.1 Scalar Field Propagator 

The quantities fo, Jo, needed to find the scalar field propagator, are seen to be 
(for a massless minimally coupled field - masses and non-minially couplings 
can be accomodated by adding £i? + m 2 to these functions) 

fo = -\z~ 2 (69) 
fo = -\z~ 2 (70) 
From which we get the function n to be equal to (remember go = c\ = fo) 
TI{X , x') = -I.^ /; ^^ir - 7 -(z- 2 + ,'-V* (71) 

As usual, the remaining coefficients r 2 ,r 3 ,... are found by differentiation of 
this quantity in appropriate ways. The remaining c n 's are found to be 

C2 = -\z~' c 3 = -\z-* (72) 

For a massless scalar field, an expression for the propagator is known already 
[|j~3|l , namely 

G '^ x ^} = A7t 2 zz , sin hf 4 (f2_ g^(t-t'Y) (73) 
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with 

2zz' 

Defining coordinates 

x = g~ 1 e g ^ sinh gr x 3 = g~ 1 e 9 ^ cosh gr 
the propagator can be rewritten as 

g\x,x) = — 9 y - 

Att 2 sinh x ^ g 2 {^ — t' + inf3) 2 — \ 2 — ie 

where (3 = ^ is an inverse temperature (appearing due to the Davies-Fulling- 
Hawking-Unruh effect, |], an d 

coshx = cosh #(£ - f) + ^ 2 a- 9(?+f,) ((xi - x[) 2 + (x 2 - x 2 ) 2 ) 

The integer n is a kind of winding number. For x 2 = x' 2 = x 3 = x' 3 = the 
geodesic distance becomes, Jl3| 

A = %- 2 sinh 2 (^(r-r')) 

For simplicity we can set x = x' = 0, y = y' = to get a function of t — t' and 
z' only (notice that neither G, nor G" are functions of z — z' but depends 
on more complicated combinations of z, z'). The quantity £4 is then simply 



£4 = Arccosh 



2 '2 
z r — z 



2zz' 

Taylor expanding around z = z' ,t = t' we get 
rv /x -1 2g 2 (t-t') 2 

( ,- 1 > ) (^ + v(«-m«-')) + 

<* - *'> 2 + 3 ^- i,)2 r,- fa + r2) ) h- ((* - « - 0") 

Since there is no singularity, this function does not satisfy the Hadamard 
condition, and can thus at most be considered as an approximation to a 
physical quantity 
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4.1.2 Gauge Boson Propagator 



In order to write down the propagator of a Yang-Mills field, we need to 
evaluate the matrices A, B, C. As we do not have a tractable expression for 
the mean fields we cannot, however, write down these explictly; only their 
curvature (or abelian) contribution can be found explicitly. For the first of 
these matrices we then get 



/ z- 



A n = - 



V 






















/ 



+ mean field terms 



(74) 



and similarly for the remaining matrices 



/ z 



B r 



























.-4 



+ mean field terms 



\ 





























+ mean field terms 



(75) 



(76) 



Thereby the coefficient A becomes 

^ 1 ^ 
' 




An 



t[ s °- > + mean field terms 



(77) 



V o o o i / 



and quite similarly for B, C 



B n = - 



nn 



9 

49 



( 1 ^ 





V o o o i / 

/ 1 \ 



V o o o i / 



_( S =0) 



+ mean field terms 



\ s °- ) + mean field terms 



(78) 



(79) 



22 



Since these matrices are all diagonal (the only possibly non-diagonal term 
comes from (gf b a c (d n A mc — d m A^ l )), i.e. one of the two terms we cannot handle 
analytically), we can actually apply the Bessel functions to them without any 
great difficulty to simply obtain the Bessel function of the (scalar) element 
multiplied by the constant 4x4 matrix appearing in A, B, C. Since the 
coefficient of proportionality between the vector and scalar case is not the 
same for all three terms (it is 3/7 for the first two, and the square of this for 
the last) we do not get just the scalar propagator multiplied by a constant 
4x4 matrix, instead we get the coefficients ot\,f3\ to be 



2(4tt)" 




3 A 

7-2 " 74^ 



x 



+mean field terms 

( 1 \ 





V 1 / 



+ mean field terms 



/ 1 \ 


V 1 / 

(80) 

(81) 



with r n given by the s = case. For an abelian field such as the Maxwell field, 
the mean field terms become irrelevant, and the expressions above becomes 
exact. 



4.2 Heat Kernel and Propagators in the Friedman- 
Robert son- Walker Metric 

An important case is the cosmologically interesting Friedman-Robertson- 
Walker space-time 

ds 2 = dt 2 - a 2 (t) (d X 2 - f (x) (d0 2 + sin 2 9d(j) 2 ) ) (82) 

in which the vierbeins can be taken to read 



1 1 ^ 








( 









( 


\ 







a 


























p a — 




af 




p a — 







V o J 




I J 




K 





I 






a f sin 9 J 



(83) 
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4.2.1 Scalar Field Propagator 

With the vierbeins given by equation (|83| ) we get 

fo = ^Qe-(9 m (ee^))) 2 -^(^ m (ee^)) 

= ^a- 2 {3a - a 2 /- 1 /' - a 2 cotO) 2 - 1(3 sin ^(a^a) - d(/'sin0 + /cos0)) 
+ l -a- 2 {-Zaf- 2 f + a 2 d x {f- 2 f)) + 1 sin (84) 
fo = ^a- 2 (?,a - a 2 f- 1 f ~ a 2 cotd) 2 + ^sined t {a- 1 a)- a{f sine + fcosd)) 
- l -a- 2 {-Zaf- 2 f + a 2 d x (f- 2 f')) - 1 sin0 (85) 

We note that fo and fo only differ in some relative signs. For the geodesic 
distance A we have not been able to find an analytical expression; in general 
one would then have to find this numericaly. 
The coefficient T\ becomes 

n = e- 2VK I e 2T \]a~ 2 (M - a 2 f- 1 f - a 2 cote) 2 + - (3sin6d 2 \na - d(f sm9 + fcos9)) 
Jo 4 v '2 



l -a- 2 {-Zar 2 f 2 + a 2 d x {r 2 f)-\ 



dr+^(Ux) + fo(x'))e- 2 ^ (86) 



For massless conformally coupled scalar fields some results are known pi 



15] in the static case; for K = — 1 one has 



^ tat '- (X ' X ' } 4vr2 C sinh( X - X ')((X ~ X') 2 + (v ~ V' ~ *t) 2 ) ^ 
while for K = +1 (the Einstein universe) one has 

G ' s ta.t,+( x > x ') = (87r 2 c(cos(r/ - V - ie) - cos(x - x'))' 1 (88) 

with i] the conformal time, and c given by i? = 6K/c. In order to compare 
these results with ours, we must first notice that in the static case with 
9 = 9', = 0' the geodesies are trivial 

t(r) = ar + b x{ T ) = ar + P 
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with these simplifications, the coefficient T\ becomes 



-2%/A 



,2t 



\al{r l f + cot0f - \d x U- 2 f) - \ sine 



dr + f e 



-2-v/A 



Now, since the geodesies are trivial, we can actually find the geodesic distance 

A = (t-t') 2 -^(x-x') 2 

where ao is some constant (the value of the scale factor). 

Taylor expanding the functions G' stat ± around rj = rj',X = x' taking 9 = 



0'. 



GL 



V by symmetry, one finds 
1 



4C7T 2 (?7 — if) 



vT, (X - X?) 



(X 



G a 



+o((v - 
-1 



X 

/\3\ 



4c7T 2 (?7-V) 4 24c7r 2 (r/-r/') 2 



/\2 



stat,- 



4cTt 2 (r] — rj') 2 

(x - X? 



48C7T 2 

1 



4c7r 2 (ry 

(x-x 



rj 

A3 



'Vl 



160C7T 2 

31(r7 - Ty') 2 " 

24cvr 2 (r7 -r]') 2 ' 2880C7T 2 ' 120960cvr 2 



+ 



11 



+ 



We notice that the singularity goes like (rj — i]')~ n , and is thus independent 
of Xi x'i we seem only to get a singularity when both rj — rj and \ — x' 
tend to zero. On the other hand, Bunch and Davies, [15[], have proven that 
the propagators obtained by a mode sum has the following structure in the 
coincidence limit x' —> x even in the non-static limit. 



a 



•i 



0(e 2 )) 



167r 2 e 2 E v 

and thus it seems to satisfy the Hadamard condition. It has been shown, 
however, by Pirk, [24], at least for spatially flat Friedman- Robertson- Walker 
space-times that adiabatic vacua are Hadamard states if and only if they are 
of infinite adiabatic order. Such subtleties are absent in our approach. 
A very general form for the propagator of a massive, conformally coupled 
scalar field in the flat Friedmann-Robertson- Walker universe has been found 



by Charach and Parker, |Tj| , and it contains the older results of Narimi and 



+ 
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Azuma, f[? |, or Chitre and Hartle, as special cases. This general formula 
is 



G(x,x') 



7T 



(2vr) 2 4tt' 



_ J e i k .(x-x0^2) (mt; 



B 



CI 



+#£V<) 



m 



terms of Hankel functions. The coefficients are restricted by 

kn — \Bk\ 2 e k7T = j, for various choices of these one then gets the 
aforementioned results by Narimi and Azuma and by Banerjee and Hartle. 
In the above formula t > is the largest of t, if , while t < is the smallest of the 
two. That this actually satisfies the Hadamard condition in this formulation 
is far from clear. Furthermore, the limit of m — > seems to be ill-defined as 
the Haenkel functions go like z~ u where z is the argument and v the order, 
contrary to this, our formula is valid for as well massive as massless fields. 
Moreover, since the expression for G(x, x') proposed in this paper does not 
involve an integral over the order of Bessel functions but only integrals of 
geometric quantities along geodesies, it might be more convenient from a 
practical view point. 

For finite temperature a result has been found by Banerjee and Mallik, [HIT]. 



4.2.2 The Gauge Boson Propagator 

In order to find the propagator of a spin 1 field we must first find the matrix 
A as defined above. Inserting the vierbeins (^) one can find 



/ 3a" 2 d 2 2a~ l f- 1 af a^dcotfl \ 

2a~ 1 f~ l af a~ 2 a 2 + 2f- 2 f' 2 f^fcotO 

rtcotfl f^f'cote a- 2 a 2 f- 2 f' 2 + cot 2 Q 

v a- 2 a 2 f- 2 f 2 + cot 2 6 ) 

+ mean field dependent terms (89) 



The matrix B is then Do = i] mn d m d n applied to this. As we do not have a 
very tractable expression for the mean fields we will only list its curvature 
contribution together with that of C in table I. 
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The coefficient A becomes similarly 



/ 3^ 



A~, 



-2y/A 



af 



a 



2^ 4 + 2^ ^cotfl 

af a 2 f 2 f f 

2 cot 9 4 cot 9 4ft + cot 2 

a / a 2 f 2 








V 





-2a/A 







+ cot 2 y 



e 2 ^r 



+ 2^n (a) + (^)) e + mean field terms 



(90) 



Again, for a static Friedman-Robertson- Walker space-time this simplifies im- 
mensely to give 



A™ = e 



1 



-2a/A 



JO 











2/" 2 /' 2 f' 1 foot 9 
f- 1 f'cot9 cot 2 9 
V 



\ 




cot 2 # J 



e 2r dr + 



(A™(x)+A™(x'))e 



-2-/A 



with (as for the scalar case) 

A = (t-ty-a 2 ( X -x) 2 

Thus the Green's function for spin one contains essentially the same integrals 
as for spin zero, albeit in a somewhat simpler configuration. 



4.3 Heat Kernel and Propagators in Schwarzschild 
Space 

The Schwarzschild metric is given by 

1 



ds 2 = h(r)dt 2 



dr 2 - r 2 dVt 



h{r) 



with 



h(r) = 1 - 



2M 



(91) 



(92) 
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where M denotes the mass of the object generating the gravitational field. 
Thus the vierbeins can be taken as 



(93) 







f °, ^ 








f \ 



























p a — 
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p a — 





V o j 




V o J 




I J 




\ r sin 9 J 



According to Chandrasekhar, [^4| , the time-like geodesies are given by 

dt 



— = E/h(r) 



dr 
dr 



dr 

2 



E 



d9 
dr 



-r- = 



dr 



L 



where E is a constant (the energy) and where L is the angular momentum. 
For the radial case L = 0, one finds 



t 



E\lr%2M 



+2M log 



-(r + sinr) + (1 - E 2 )t 

tan \t h + tan |r 
tan \t h — tan \r 



where 



r cos — T 



2M 



1 -E 2 

is the original position (the geodesies describe matter falling into the black 
hole, or whatever is creating the gravitational field) and where 

th = 2Arcsin E 

Null geodesies, describing massless particles, are somewhat easier. The radial 
solution is 



t 
r 



±n(r) + c± 
±Et + c± 
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with 

U =r + 2Mlog (^J-l 
and where c± are constants. 

4.3.1 Scalar Field Propagator 

For the scalar field we only need to evaluate the two quantities /o, /o- With 
the vierbeins as given above, this is an easy task and we get 

f = r~ 3 {r -2M)+ ] -M 2 r~ 3 {r -2M)" 1 + Mr~ 3 + -r~ 2 

4r Z 

-Mr~ 4 (r - 2M)- 1 + ^M 2 r" 4 (r - 2M)~ 2 + M(r - 2M)' 1 (94) 

f = r- 3 {r -2M) + ^M 2 r*{r -2M)- 1 + Mr' 3 - -r~ 2 

4 2 

+Mr" 4 (r - 2M)~ l - ^M 2 r~\r - 2M)~ 2 - M(r - 2M)" 1 (95) 

2 



which can then be inserted into the definition of T\ 

n = e- 2 ^ C foe^dr + ±(/ (r) + / (r'))e 
Jo 2 



for a radial null-geodesic, the integral over r can actually be carried out and 
we get 



/ U{r{r))e 2T dT = (cE)~ 2 (c(E — 2M) — EM) + (E(c + E))~ 2 e 2 (E 2 + c(E — 2M) — 3EM) + 
Jo 

2E- 3 e- 2cE ~\E - 2M)(Ei(2 + 2cE~ x ) - Ei(2cE~ 1 )) + ^M 2 x 

-c{E - c - 2M) - -E(E - 2M) - 2M 2 + -e 2 c(c + E + 2M) + 
2 4 2 

^-e 2 £(£ + 2M) + 2M 2 e 2 + 8E' 1 M 3 e' 2{c ~ 2M)/E x 



(Ei(2(c + E- 2M) /E) - Ei(2(c - 2M) / E))] 
M 

(Ei(2c/£) - Ei(2 + 2c/£))] 



-c-'E- 2 - -e 2 E~ 2 (2c + 3E)(c + E)~ 2 + 2E- 3 e~ 2c ' E x 

Z Z 
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''cEY 1 - e 2 E-\c + E)~ l + 2e- 2c/£ E- 2 (Ei(2 + 2c/ E) - Ei(2c/£)) 



+ 



'" [3c 2 E{E + 2M) + 3cE 2 M + 4M 2 {E 2 + 2c 2 + cE) 



2 

M 

L 24c 3 £ 3 M 3 

e 2 (24£ 3 M 3 (c + E) 3 ) -1 (3c 2 £(£ + 2£ 2 + 2M) + 3£ 4 + E 3 M(9 + 15c) + 

4M 2 (4£ 2 + 2c 2 + 5c£?)) + (12M 3 E 4 )- 1 e - 2c/B (3 J B 2 + 6£M + 8M 2 ) x 

(Ei(2c/£) - Ei(2 + 2c/£)) + 

(\QEM A )' l e~ 2c/E (Ei(2c/£) - Ei(2 + 2c/E) + 

Ei(2(c + E-2M) /E)e~ 4M/E - Ei(2(c - 2M) /E)e- 4M/£ )" - 

-M 2 x [(12EM 2 c 3 )^ - (16£(2M - c)M 4 ) _1 + 

e 2 (f 6£(2M — E — c)M 4 ) -1 - e 2 (12EM 2 (c + E) 3 ) -1 + 

3(SM 5 ) _1 e" 2c/B (Ei(2 + 2c/£) - Ei(2c/£)) + 
8 

(48c 2 E 3 M 4 )" 1 (c(9£ 2 + 12EM + 8M 2 ) + 6E 2 M + AEM 2 ) + 
e 2 (48£ 3 (c + EfM 4 )- 1 (9E 2 (c + E + 2M) + \2cEM + 4M 2 (3£ + 2c)) + 
(48 J B 4 M 4 )- 1 e- 2c / £ (Ei(2c/£) - Ei(2 + 2c/ E)) + 

(8E 2 M 5 )- 1 e - 2 ^- 2M ^ E (Ei(2(c + E-2M) / E) - Ei(2(c - 2M) / E)) - 
ME -i e -2(c-2M)/E ( Ei ( 2 (c + E-2M) /E) - Ei(2(c - 2M) /E)\ 

which looks very complicated but is actually just a combination of polyno- 
mials, exponentials and exponential integral functions. 
Unfortunately, the geodesic distance A is rather difficult to find analytically 
in Schwarzschild space-time, only radial geodesies are easily obtained. In 
many aplications, though, these actually suffice. Along these geodesies the 
differential operators appearing in the expressions for t 2 ,t 3 are 

'cV [>l r) dr 



d P = e;x^ = (h~ i m-h 1 / 2 r)^- (96) 



□o = V pq d p d g = (h- 1 ii-hrr)— + (h~ 1 i 2 -hr 2 )— (97) 

or or 

Fairly simple expressions. But we will not write down the explicit form for 
Tii r 3 which follows from the application of these operators as the results are 
too involved. 
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Candelas and Jensen, SO] , have also found a propagator for the Schwarzschild 
background. They use the Hartle-Hawking vacuum (i.e. the Green's func- 
tions is demanded to vanish as the spatial distance tends to infinity and to 
be periodic in imaginary time), with this they arrive at 



G(—ir, r, 9, 0, —it', r', 9', 0') 



32tt 2 M 2 



£(2Z + l)P i (cos 7 )P^<)Q/(£>) + 



l=o 



OO 1 oo 

£ _ C0S (nK(r - /)) 52(21 + l)P l (cos -y)pttt<)q?& 

71=1 n 1=0 



where Pi, Qi are Legendre polynomials, pf, qf are solutions to 

f = o 



at; 



where £ = j? — 1 and 



cos 7 



16(£ 2 - 1) 



max(£, f ') 

cos 9 cos 0' + sin 9 sin cos( 



They prove that 

lim G(x, x') 



8ir 2 A(x,x') ' 12(8vrM) 2 l-^f ' (8ttM) 



+ 



with F(r) small. Like our results, then, Candelas and Jensen finds a Green's 
function with only the A _1 -singularity, whether they also get a A log A term 
is, however, not so clear, since not much is known about q^pf and the sums 
involving these and/or Legendre polynomials is difficult to actually carry out. 
If one calculate the finite part of our propagator as x' — *> x, one finds 

M 



32tt 2 



(48M 3 - M 2 r (48 + 52M 2 ) + Mr 2 (19 + 159M 2 ) - 158M 2 r 3 + 
SMr 4 - r 5 (10 - 4M 2 ) - 2Mr 6 ) (8r 7 (r - 2M) 3 )" 1 + 
— (48M 3 - M 2 r(48 + 52M 2 ) + Mr 2 (19 + 159M 2 ) - 158M 2 r 3 + 
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• 5 (10 - 4M 2 ) - 2Mr 6 ) (24r 10 (r - 2M) 6 )" 1 - 
X — (240M 4 - M 3 r(384 + 216M 2 ) + M 2 r 2 (198 + 666M 2 ) - Mr 3 (30 + 793M 2 ) + 
461MV - 132Mr 5 + r 6 (15 - 4M 2 ) + 2Mr 7 ) (8r 7 (r - 2M) 3 )" 1 



66Mr 4 
1 

32^ 



i.e. one finds more terms which diverge as one approaches the Schwarzschild 
radius. A discrepancy likely to be due to (1) the different vacua (we do not 
have an explicit temperature) and (2) the different orders of the respective 
approximations. 



4.3.2 Gauge Boson Propagator 

The curvature contribution to the matrix A is in this case, with the vierbeins 
as given above, easy to find and we get 



A n — 





/ h^drh- 1 / 2 ) 2 
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1 





h 2 {d r h- 1 / 2 ) 2 + 2 A 


h 1 / 2 cote 





2 





h 1 / 2 cote 


h+cot e 







V o 








fc+cot 9 
r 2 / 



+ mean field terms 
whereby the matrices B, C become as follows 



(98) 



B n = 



/ 2M 2 (27M 2 -32Af r+10r 2 ) 
r 6 (r-2M) 2 









V 



o 





2(243Af 4 -400M 3 r+240M 2 r 2 -62Mr 3 +6r 4 
r 6 (r-2M) 2 

(2M-r)(10M-4r-10M cos 20+3r cos 20) cot 9 esc 2 t 





(2M-r)(10M-4r-10M cos 20+3r cos 20) cot 6 esc 2 9 
54M 2 -38Afr+6r 2 _ (14M-6r) cot 29 , cot 9 esc 2 9 



+ mean field terms 







2r 2 -2M 2 -2Mr _ (14M-6r) cot 9 . cot esc 2 9 
„2 „5 ~T „4 



(99) 



where we, for reasons of space, have split-up the matrix B and 
1 



12 



x 
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n 4(27M 3 -34M 2 r+14Mr 2 -2r 3 ) 2 (5M-2r) 2 cot 2 6> , (r-2M) esc 4 9 n 

U r 7 (r-2M) 5 r 6 (r-2M) 2 + r 7 U 

n (5M-2r) 2 cot 2 , (r-2M)csc 4 6> 6M-2r _ 4cot 2 6> , esc 4 6> r, 

U r-6(r-2Af) 2 r 7 r 4 r 5 (r-2M) r 6 U 

r, r, n 6Af-2r 4 cot 2 g , esc 4 g 

U U U r 4 r 5( r _2M) "T" r 8 / 

+ mean field terms (100) 
And this can then be inserted into the formula fl3Tj| ) for the propagator 



4.4 Conformally Flat Space-Times 

As our last example we will consider a conformally flat space-time with line 
element 

ds 2 = C{r)f{drf - dijd^dx 1 ) (101) 

This could for instance be de Sitter space or a general spatially flat Friedman- 
Robert son- Walker metric in arbitrary dimensions. For Friedman-Robertson- 
Walker space-time the conformal time 77 is defined to be /* a~ x dt and C(rj) 
is just the usual scale factor a(t) expressed in this new coordinate, whereas 
for de Sitter space-time 77 = — aexp(— t/a) where a = J\2/R with R the 
curvature scalar (i.e. a constant in this case). The calculations will be carried 
out in d — 4 dimensions, but can of course easily be generalized to any d. 
For the metric as given above, we can choose the vierbeins to be 
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f ^ 




( \ 



























p a — 

e 2 _ 


c(v) 


p a — 





V J 




V j 




V J 




V C(rj) J 



(102) 

4.4.1 Scalar Boson Case 

The d'Alembertian can be written as 

□ = D + S a d a 



(103) 
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where C = d^C. We then have 

fo = 2 -lc 2 C-*- 3 -CC~ 3 (104) 

; = ^c 2 c- 4 + ^cc- 3 (105) 

whereby we get 

n = e" 2 ^ jf 1 (^(7 2 ^ 4 + ^^" 3 ) e 2T dr 

+\ {^C^fC{vr A + + 3C( V )C( V )- 3 + 3C(t/)C(V)- 8 ) e- 2 ^ 

(106) 

and the coefficients C2, C3 are found to be 

c 2 = -^(C'- 5 C' (4) + C- 6 (C' + l)C^)-3C- 6 (7 2 + ic- 7 (51C' + 174)C Y2 C 

-21C- 7 C 4 - ^C~ 8 C 4 (107) 

147 . .. ... 

c 3 = ^C" 10 C 2 C 2 - 147C- n C 4 C + 147C- 12 C- 6 (108) 

It turns out that in the particular case of de Sitter space-time, another for- 
mula for the propagator of a scalar field is known in the literature |T], [2"2] ], 
namel 



r'( >\ 1 ( l ^ v( 3 ^ 3 9 1 i (r/-r/-26) 2 -(x-x')% 

G x, x = — ir) sec7rz/F — hzA — v\ 2; H 

v ' ; 167ra 2V 4 ; v 2 '2 ' ' 4777/ ; 

(109) 

with 

v 2 = - Uim'R- 1 + 

Clearly this can not be the final answer as it diverges for a massless minimally 
coupled field (in this case v = |). Let us nonetheless compare it with our 
result obtained by the method described in this paper. The first thing we 



8 For thermal fields, a propagator was also found by Chardury and Mallik, p3[ , but it 
will not be discussed here. 
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notice is that, even on the level of series expansion as in equation fl2"2|), our 
expression is manifestly divergent as x — > x' (due to the A -1 , log A terms), 
whereas the expression ( |109| ) only diverges when all the terms in the series 
expansion of the hypergeometric function is included (it then goes to T(— 1)). 
More explicitly, by letting x — > x' along a non-null geodesic with tangent t M , 
tH^ = £ = ±1, one has (for d = 2, but the form is the same in any 
dimension, only the coefficients change) 



G(x, x') 



-(4?r 



2 7 + In ie 2 Si? + V(~ +fi) + - /i) + 



-e 2 T,R - e 2 Z{£R + m 2 ) (27 + In ~e 2 ZR+ 



with 



2£ — m 2 a 2 



Now, in this case A = |Se 2 , so what we capture here is only the log A diver- 
gence and not the more common A -1 one, consequently it does not satisfy 
the Hadamard condition. The result ( |109|) follows from a mode summation 
using an adiabatic vacuum, the failure to meet the Hadamard condition is 
then likely to be a result of not going to infinite adiabatic order. The adia- 
batic vacuum has the advantage of being very similar to flat space-time, but 
is, alas, only possible in some space-times. The Green's function following 
from the heat kernel as put forward here, on the other hand, comes with 
a "natural vacuum", which apparently works in all cases. The difference 
between our result and (|109|) can thus come from many sources: 



Ambiguity in the definition of a vacuum 
which does not always work). 



uses adiabatic vacuum 



Ambiguity in the definition of a propagator (the Green's function as 
defined in this paper is rather unique, and from this all propagators 
- positive frequency, advanced, retarded, Whightmann etc. - can be 
found) . 

The approximations made (our result is highly non-pertubative and 
seems to include more terms than in ( |109| ). 
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This also summarizes the differences found in the other specific geometries. 



4.4.2 Gauge Boson Case 

The tensor 

g r = rT^e* - d q e» n )el 
is found to have the following non-vanishing components only 

£« = Ef = -CC~ X (110) 

with no summation over the index i = 1, 2, 3 implied. From this we get the 
mean-field independent part of A to be 

= -^C 2 C- 2 + mean field terms (111) 



whereby 



~Jq 6 ™ 6 ~ 2 ^ (pWfClri)- 2 + C{rf) 2 C(r}')- 2 ) + mean field(sll2) 
For the matrix-coefficients B, C we get similarly 

B™ = -5™ ( 2C~ 4 C 2 + 2C~ 4 ct - 10C- 5 C 2 C + 6C~ 6 C 4 + 2C 2 CC~ 4 - 2C~ 5 C 4 



+ mean field terms (113) 
C™ = — 5™ (CCC- 2 -C- ?> C ? ') 2 + mean field terms (114) 

And once more one can get the remaining terms by simple differentiation of 
these. 



5 Conclusion and Outlook 

We have developed a way to determine the propagators of quantum fields 
in curved space-time explicitly. It is, for practical reasons, presently only 
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possible to calculate an approximation to the propagator. This approxima- 
tion, while excellent for scalar and spinor fields and relatively good for spin 1 
gauge bosons is, in the latter case, particularly for non-abelian gauge fields, 
rather complicated to work with so one might hope to change this by altering 
the approach by which we found this expression. In all the cases, though, 
the practical calculation of the Green's function, was reduced to that of cal- 
culating a few functions (for spin zero and one half, f , / ) or three matrices 
(for vector bosons, A, B, C) upto interactions with other fields than the grav- 
itational background. 

It turned out, however, that our non-perturbative propagators, though sat- 
isfying the Hadamard condition and having a very simple physical meaning, 
differed from what other authors had found in various cases. This differ- 
ence can be traced back to the uniqueness of the vacuum as presented here, 
and which differs, in general, from the adiabatic or conformal vacuum usu- 
aly used. The ambiguity in the definition of a vacuum is related to the 
insufficiency of the naive particle concept in curved space-times. By always 
referring to a local comoving observer, we seem to have arrived at a better 
defined particle concept, since we are then able to relate to a flat space-time 
at each point along the trajectory, and of course, in flat space-time there 
is no problem with the particle concept, as the general success of ordinary 
quantum field theory shows. It should furthermore be emphasized that the 
method put forward in this paper is very general, as we do not to assume any 
particular feature of the metric. If the space-time under consideration does 
show some particular features (such as weak or slowly varying curvature, 
asymptotic or conformal flatness and so on), then of course, the calculation 
of the propagator by our method can be simplified, or rather adjusted, by 
then only including the necessary terms in the general expansion of the heat 
kernel. As these coefficients are actually rather easy to find this adjustment 
is straightforward in each particular case. 

Looking at our results for, say, Friedman- Robert son- Walker and Schwarzschild 
geometries, we notice the appearance of terms which survive if the flat space 
limit is understood as a(t) — > oo and M — > respectively. These effect 
are most likely due to the difference in topologies; for Friedman-Robertson- 
Walker space-times we still have K = ±1,0 even in the limit a — > oo, and 
in the Schwarzschild case, the limit M — > corresponds to Minkowski space 
with the origin removed. It is not at present clear what the physical impli- 
cations of such corrections are, however. When doing perturbation theory 



37 



in flat space-time, one should of course in principle take into account that 
we do not live in Minkowski space, and thus use, say, the a — > oo limit of a 
Friedman-Robertson- Walker propagator. One would expect any such effect 
to be very small, but it is not a priori certain as the results presented here 
are in deed non-perturbative, and even though the perturbative effect is very 
small, the true, non-perturbative effect need not be. Given the present day 
limitations on accelerator power, the only chance of seeing a quantum grav- 
itational effect, or an effect coming from the coupling of the quantum fields 
to the space-time geometry, is either by studying the early universe and/or 
black holes or to look for some non-perturbative effect. 
Naturally, the method of calculating the heat kernel as presented here in this 
paper is in no way restricted to quantum field theory in curved space-time. 
It should also be of considerable mathematical interest, as it is known that 
the heat kernel of the Laplace-Beltrami operator contains information about 
the large scale structure of the underlying manifold. Non-perturbative quan- 
tum field theory in ordinary Minkowski space is of course just a particular 
example, and the uses there include quark-gluon plasma. 
The uses for the propagators include refinements in determining the spectrum 
of Hawking emission from black holes as well as the emission of gravitational 
waves from the space-time singularity (the Big Bang), applications that we 
hope to entertain in subsequent papers. Another application, which we have 
not yet considered, is the evaluation of the trace-anomaly of the energy- 
momentum tensor in general. Lastly we should note that, having obtained 
expressions for the propagators in this paper and for the one-loop effective 
actions in a previous one, we are now in a position to find the effective action 
to any loop order (at least formally). A paper on this is in progress. 
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Table 1: The mean field independent contribution to the symmetric matrices 
B and C for a Friedman-Robertson- Walker space-time. 
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